% function [baseMVA, bus, gen, branch, genpara, NTswitch, sbspara,Hess_coef,u0,mu_coef] = case30
baseMVA = 100.00;

% bus  type1 type2  Pdemand  Qdemand   Gshunt   Bshunt   Vmodule    Vangle  baseKV  Vmax  Vmin
% type1: PQ bus = 1; PVbus = 2; reference bus = 3; isolated bus = 4
% type2: non-load = -1; load -- no shedding = 0; load -- shedding =1
bus = [
     1     3    -1      0.0000      0.0000      0.0000      0.0000      1.0000      0.0000    135.0000      1.0500      0.9500
     2     2     1     21.7000     12.7000      0.0000      0.0000      1.0000     -0.4150    135.0000      1.1000      0.9500
     3     1     1      2.4000      1.2000      0.0000      0.0000      0.9831     -1.5220    135.0000      1.0500      0.9500
     4     1     1      7.6000      1.6000      0.0000      0.0000      0.9801     -1.7950    135.0000      1.0500      0.9500
     5     1    -1      0.0000      0.0000      0.0000      0.1900      0.9824     -1.8640    135.0000      1.0500      0.9500
     6     1    -1      0.0000      0.0000      0.0000      0.0000      0.9732     -2.2670    135.0000      1.0500      0.9500
     7     1     1     22.8000     10.9000      0.0000      0.0000      0.9674     -2.6520    135.0000      1.0500      0.9500
     8     1     1     30.0000     30.0000      0.0000      0.0000      0.9606     -2.7260    135.0000      1.0500      0.9500
     9     1    -1      0.0000      0.0000      0.0000      0.0000      0.9805     -2.9970    135.0000      1.0500      0.9500
    10     1     1      5.8000      2.0000      0.0000      0.0000      0.9844     -3.3750    135.0000      1.0500      0.9500
    11     1    -1      0.0000      0.0000      0.0000      0.0000      0.9805     -2.9970    135.0000      1.0500      0.9500
    12     1     1     11.2000      7.5000      0.0000      0.0000      0.9855     -1.5370    135.0000      1.0500      0.9500
    13     2    -1      0.0000      0.0000      0.0000      0.0000      1.0000      1.4760    135.0000      1.1000      0.9500
    14     1     1      6.2000      1.6000      0.0000      0.0000      0.9767     -2.3080    135.0000      1.0500      0.9500
    15     1     1      8.2000      2.5000      0.0000      0.0000      0.9802     -2.3120    135.0000      1.0500      0.9500
    16     1     1      3.5000      1.8000      0.0000      0.0000      0.9774     -2.6440    135.0000      1.0500      0.9500
    17     1     1      9.0000      5.8000      0.0000      0.0000      0.9769     -3.3920    135.0000      1.0500      0.9500
    18     1     1      3.2000      0.9000      0.0000      0.0000      0.9684     -3.4780    135.0000      1.0500      0.9500
    19     1     1      9.5000      3.4000      0.0000      0.0000      0.9653     -3.9580    135.0000      1.0500      0.9500
    20     1     1      2.2000      0.7000      0.0000      0.0000      0.9692     -3.8710    135.0000      1.0500      0.9500
    21     1     1     17.5000     11.2000      0.0000      0.0000      0.9934     -3.4880    135.0000      1.0500      0.9500
    22     2    -1      0.0000      0.0000      0.0000      0.0000      1.0000     -3.3930    135.0000      1.1000      0.9500
    23     2     1      3.2000      1.6000      0.0000      0.0000      1.0000     -1.5890    135.0000      1.1000      0.9500
    24     1     1      8.7000      6.7000      0.0000      0.0400      0.9886     -2.6310    135.0000      1.0500      0.9500
    25     1    -1      0.0000      0.0000      0.0000      0.0000      0.9902     -1.6900    135.0000      1.0500      0.9500
    26     1     1      3.5000      2.3000      0.0000      0.0000      0.9722     -2.1390    135.0000      1.0500      0.9500
    27     2    -1      0.0000      0.0000      0.0000      0.0000      1.0000     -0.8280    135.0000      1.1000      0.9500
    28     1    -1      0.0000      0.0000      0.0000      0.0000      0.9747     -2.2660    135.0000      1.0500      0.9500
    29     1     1      2.4000      0.9000      0.0000      0.0000      0.9796     -2.1280    135.0000      1.0500      0.9500
    30     1     1     10.6000      1.9000      0.0000      0.0000      0.9679     -3.0420    135.0000      1.0500      0.9500
];

% bus       Pg         Qg     Vsetpoint  base   status  Pmax  Pmin  Qmax  Qmin
gen = [
     1     25.9738     -0.9985      1.0000    100.0000     1     1     80.0000      0.0000    150.0000    -20.0000
     2     60.9700     31.9990      1.0000    100.0000     1     1     80.0000      0.0000     60.0000    -20.0000
    22     21.5900     39.5700      1.0000    100.0000     1     1     50.0000      0.0000     62.5000    -15.0000
    27     26.9100     10.5405      1.0000    100.0000     1     1     55.0000      0.0000     48.7000    -15.0000
    23     19.2000      7.9510      1.0000    100.0000     1     1     30.0000      0.0000     40.0000    -10.0000
    13     37.0000     11.3529      1.0000    100.0000     1     1     40.0000      0.0000     44.7000    -15.0000
];

% fbus  tbus    r        x        b      ratea    rateb     ratec Vf/Vt-ratio  angle  status
% rateA: long term rating; rateB: short term rating; rateC: emergency rating
branch = [
     1     2      0.0200      0.0600      0.0300    130.0000    130.0000    130.0000      0.0000      0.0000     1
     1     3      0.0500      0.1900      0.0200    130.0000    130.0000    130.0000      0.0000      0.0000     1
     2     4      0.0600      0.1700      0.0200     65.0000     65.0000     65.0000      0.0000      0.0000     1
     3     4      0.0100      0.0400      0.0000    130.0000    130.0000    130.0000      0.0000      0.0000     1
     2     5      0.0500      0.2000      0.0200    130.0000    130.0000    130.0000      0.0000      0.0000     1
     2     6      0.0600      0.1800      0.0200     65.0000     65.0000     65.0000      0.0000      0.0000     1
     4     6      0.0100      0.0400      0.0000     90.0000     90.0000     90.0000      0.0000      0.0000     1
     5     7      0.0500      0.1200      0.0100     70.0000     70.0000     70.0000      0.0000      0.0000     1
     6     7      0.0300      0.0800      0.0100    130.0000    130.0000    130.0000      0.0000      0.0000     1
     6     8      0.0100      0.0400      0.0000     32.0000     32.0000     32.0000      0.0000      0.0000     1
     6     9      0.0000      0.2100      0.0000     65.0000     65.0000     65.0000      0.0000      0.0000     1
     6    10      0.0000      0.5600      0.0000     32.0000     32.0000     32.0000      0.0000      0.0000     1
     9    11      0.0000      0.2100      0.0000     65.0000     65.0000     65.0000      0.0000      0.0000     1
     9    10      0.0000      0.1100      0.0000     65.0000     65.0000     65.0000      0.0000      0.0000     1
     4    12      0.0000      0.2600      0.0000     65.0000     65.0000     65.0000      0.0000      0.0000     1
    12    13      0.0000      0.1400      0.0000     65.0000     65.0000     65.0000      0.0000      0.0000     1
    12    14      0.1200      0.2600      0.0000     32.0000     32.0000     32.0000      0.0000      0.0000     1
    12    15      0.0700      0.1300      0.0000     32.0000     32.0000     32.0000      0.0000      0.0000     1
    12    16      0.0900      0.2000      0.0000     32.0000     32.0000     32.0000      0.0000      0.0000     1
    14    15      0.2200      0.2000      0.0000     16.0000     16.0000     16.0000      0.0000      0.0000     1
    16    17      0.0800      0.1900      0.0000     16.0000     16.0000     16.0000      0.0000      0.0000     1
    15    18      0.1100      0.2200      0.0000     16.0000     16.0000     16.0000      0.0000      0.0000     1
    18    19      0.0600      0.1300      0.0000     16.0000     16.0000     16.0000      0.0000      0.0000     1
    19    20      0.0300      0.0700      0.0000     32.0000     32.0000     32.0000      0.0000      0.0000     1
    10    20      0.0900      0.2100      0.0000     32.0000     32.0000     32.0000      0.0000      0.0000     1
    10    17      0.0300      0.0800      0.0000     32.0000     32.0000     32.0000      0.0000      0.0000     1
    10    21      0.0300      0.0700      0.0000     32.0000     32.0000     32.0000      0.0000      0.0000     1
    10    22      0.0700      0.1500      0.0000     32.0000     32.0000     32.0000      0.0000      0.0000     1
    21    22      0.0100      0.0200      0.0000     32.0000     32.0000     32.0000      0.0000      0.0000     1
    15    23      0.1000      0.2000      0.0000     16.0000     16.0000     16.0000      0.0000      0.0000     1
    22    24      0.1200      0.1800      0.0000     16.0000     16.0000     16.0000      0.0000      0.0000     1
    23    24      0.1300      0.2700      0.0000     16.0000     16.0000     16.0000      0.0000      0.0000     1
    24    25      0.1900      0.3300      0.0000     16.0000     16.0000     16.0000      0.0000      0.0000     1
    25    26      0.2500      0.3800      0.0000     16.0000     16.0000     16.0000      0.0000      0.0000     1
    25    27      0.1100      0.2100      0.0000     16.0000     16.0000     16.0000      0.0000      0.0000     1
    28    27      0.0000      0.4000      0.0000     65.0000     65.0000     65.0000      0.0000      0.0000     1
    27    29      0.2200      0.4200      0.0000     16.0000     16.0000     16.0000      0.0000      0.0000     1
    27    30      0.3200      0.6000      0.0000     16.0000     16.0000     16.0000      0.0000      0.0000     1
    29    30      0.2400      0.4500      0.0000     16.0000     16.0000     16.0000      0.0000      0.0000     1
     8    28      0.0600      0.2000      0.0200     32.0000     32.0000     32.0000      0.0000      0.0000     1
     6    28      0.0200      0.0600      0.0100     32.0000     32.0000     32.0000      0.0000      0.0000     1
];

% bus   model  startup  shutdown n_coef   Xd'       H=1/2*Tj  D
% model: 1 -- piecewise linear   2 -- polynomial
genpara = [
     1     2      0.0000      0.0000     3      1.2100      3.0000      0.0000
     2     2      0.0000      0.0000     3      1.5500      2.0000      0.0000
    22     2      0.0000      0.0000     3      1.3750      1.2500      0.0000
    27     2      0.0000      0.0000     3      1.3125      4.3750      0.0000
    23     2      0.0000      0.0000     3      1.5250      8.7500      0.0000
    13     2      0.0000      0.0000     3      1.5000      5.0000      0.0000
];

% contingency set
% Type    Loctn1   Loctn2    Occurr      Clear_fault    Clear_load     R           X
NTswitch = [
  1110     2     5      0.0000      0.1800      0.3000  1.0000e-010  1.0000e-010
];

% parameters for transient stability
sbspara = [
        0.01
        110 *pi/180
       60.00
];
